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(57) ABSTRACT 


Real time cardiac electrical data are received from a patient, 
manipulated to determine various useful aspects of the ECG 
signal, and displayed in real time in a useful form on a com- 
puter screen or monitor. The monitor displays the high fre- 
quency data from the QRS complex in units of microvolts, 
juxtaposed with a display of conventional ECG data in units 
of millivolts or microvolts. The high frequency data are ana- 
lyzed for their root mean square (RMS) voltage values and the 
discrete RMS values and related parameters are displayed in 
real time. The high frequency data from the QRS complex are 
analyzed with imbedded algorithms to determine the pres- 
ence or absence of reduced amplitude zones, referred to 
herein as “RAZs”. RAZs are displayed as “go, no-go” signals 
on the computer monitor. The RMS and related values of the 
high frequency components are displayed as time varying 
signals, and the presence or absence of RAZs may be simi- 
larly displayed over time. 

9 Claims, 16 Drawing Sheets 
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REAL-TIME, HIGH FREQUENCY QRS 
ELECTROCARDIOGRAPH WITH REDUCED 
AMPLITUDE ZONE DETECTION 

This application is a continuation of prior application Ser. 5 
No. 09/906,013, filed Jul. 12, 2001, nowU.S. Pat No. 7,113, 
820. 

The invention described herein was made by employees of 
the United States Government and may be manufactured and 
used by or for the government for government purposes with- 10 
out payment of any royalties thereon or therefore. 

FIELD OF THE INVENTION 

The present invention relates generally to the field of elec- 1 5 
trocardiography, and more particularly to a real-time process- 
ing system and method to analyze and display electrocardio- 
graphic signals. 

BACKGROUND OF THE INVENTION 20 

Diagnosis of abnormal cardiac conditions has relied in the 
past on visible alterations in the P, QRS, and T-waves, i.e. 
portions of the electrocardiograph periodic signal. The elec- 
trocardiograph signal includes a low frequency portion and an 25 
impressed or imbedded high frequency portion, and it has 
been found that although the higher frequency portion of the 
signal is not particularly visible, it contains information that 
provides greater sensitivity in determining certain abnormali- 
ties, notably myocardial ischemia and infarction. 30 

The conventional electrocardiogram (ECG) can be a very 
insensitive diagnostic tool. For example, a significant per- 
centage of individuals presenting to a hospital emergency 
room with an actual myocardial infarction (heart attack) will 
have a normal 12 -lead conventional ECG. In addition, the 35 
conventional ECG accurately reflects only the predominant 
low-frequency electrical activity of the heart. It tells the cli- 
nician little or nothing about the less predominant high fre- 
quency components of the heart’s electrical signal embedded 
within the various lower- frequency waves of the conventional 40 
ECG. 

From off-line studies, it is known that a diminution of the 
higher frequency components within the central portion of 
the QRS complex of the ECG can be a highly sensitive indi- 
cator for the presence of myocardial ischemia or infarction, 45 
more sensitive, for example, than changes in the ST segment 
of the conventional low-frequency electrocardiogram. How- 
ever, until now, there has been no device capable of display- 
ing, in real time, changes in these high frequency QRS com- 
ponents in the monitored patient. While academic software 50 
programs have been designed that analyze the central high 
frequency QRS components, all such programs involve labo- 
rious off-line calculations and post-processing, and therefore 
have little if any clinical utility, being strictly research tools. 
Thus, there remains a need for a system and method that 55 
analyzes high frequency components over the entire QRS 
interval in real time for usefulness in the clinical environment. 
Such a system should perform, in real time, all of the complex 
digital sampling, averaging, and filtering that is required to 
generate high frequency QRS ECG signals. The system 60 
should also thereafter update these high frequency QRS ECG 
signals, as well as other derived parameters, in real time on a 
beat-to-beat basis, supplementing the diagnostic information 
being obtained from the conventional (i.e. low frequency) 
ECG complexes at the same time. 65 

The higher frequency signals in the central portion of the 
QRS ECG complex that have generated the most research 
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interest in terms of off-line detection of ischemia and infarc- 
tion are those signals in the range of 1 50 to 250 Hz. The raw, 
analog ECG signal is typically sampled at ^500 samples per 
second (to digitize the signal) in order to adequately satisfy 
the Nyquist rate of sampling at least twice the highest fre- 
quency of interest and in order to retain the information in the 
signal without loss. In the past, the sampled data have been 
stored, and then later processed to provide potentially uselul 
information to the researcher. 

On the other hand, Simpson, in U.S. Pat. No. 4,422,459, 
teaches a system which analyzes only the late portion of the 
QRS interval and early portion of the ST segment, and in an 
off-line fashion (i.e. from previously stored data) to indicate 
cardiac abnormalities, in particular the propensity for cardiac 
arrhythmia. The late portion of a post myocardial infarct 
patient’s QRS waveform contains a high frequency (40-250 
Hz) signal tail which is indicative of a tendency toward ven- 
tricular tachycardia. The system in Simpson digitally pro- 
cesses and filters a patient’s QRS signals in a reverse time 
manner to isolate the high frequency tail and avoid the filter 
ringing which would otherwise hide the signal. Thus, in order 
to do so, Simpson presupposes that the data are stored so that 
they can be processed in reverse time order. 

Albert et al. U.S. Pat. No. 5,117,833, partially focuses on 
analyzing signals within the mid-portion of the QRS interval 
for the indication of cardiac abnormality. The system of 
Albert et al. uses a known technique of building up data points 
to derive an average of heartbeat characteristics in order to 
enhance signal to noise ratio. Data are collected and filtered 
and then stored for subsequent analysis. Thus, the system 
does not teach a cardiac monitor which provides the data 
analysis immediately from the data derived from a patient, i.e. 
in “real-time”. 

Albert et al., U.S. Pat. No. 5,046,504, similarly teaches the 
acquisition of QRS data and subsequent analysis. Routine 
calculations are performed from the data previously calcu- 
lated and stored. Further, this system teaches producing a set 
of digital spectrum values representative of an approximate 
power density spectrum at each of a large number of generally 
equally spaced sampling time intervals of the ECG wave- 
form. 

Seegobin, in U.S. Pat. Nos. 5,655,540 and 5,954,664, pro- 
vides a method for identifying coronary artery disease. The 
method relies on a database of high and low frequency ECG 
data taken from known healthy and diseased subjects. Com- 
parison of the data has led to a “Score” component, indicating 
deviation of a patient’s data from the norm. This reference is 
rather calculation intensive, and does not suggest monitoring 
the condition of a patient, but rather is utilized as an off-line 
diagnostic tool. 

Hutson, U.S. Pat. No. 5,348,020, teaches a technique of 
near real-time analysis and display. The technique includes 
inputting ECG data from multiple, sequential time intervals 
and formatting those data into a two-dimensional matrix. The 
matrix is then decomposed to obtain corresponding singular 
values and vectors for data compression. The compressed 
form of the matrix is analyzed and filtered to identify and 
enhance ECG signal components of interest. As with other 
systems, this reference focuses on late potentials, a fraction of 
the QRS interval, as the tool to identify cardiac disease. 

Finally, High-Frequency Electrocardiogram Analysis of 
the Entire QRS in the Diagnosis and Assessment of Coronary 
Artery Disease by Abboud ( Progress in Cardiovascular Dis- 
eases, Vol. XXXV, No. 5 (March/April), 1993: pp 311-328) 
teaches the concept of “reduced amplitude zone” (RAZ) as a 



US 7,539,535 B1 


3 

diagnostic tool. However, this reference also uses post-pro- 
cessing, and provides no teaching of a real-time analysis 
system. 

Thus, there remains a need for an electrocardiograph that 
analyzes, in real time, the high frequency components of the 5 
QRS complex in order to provide an effective monitor for 
patients with specific cardiac function abnormalities. The 
present invention is directed to such an electrocardiograph. 

SUMMARY OF THE INVENTION 1° 

The present invention addresses these and other needs in 
the art by providing a real time display of various aspects of 
the QRS complex. The invention also provides a system for 
such a display, and a method of displaying such aspects. 

The present invention advances the state of the art by taking 
high frequency electrocardiographic data immediately as 
they are sensed from a patient, manipulating the data in con- 
junction with the conventional low frequency signals, and 2Q 
displaying the high frequency data in real time in a useful 
form on a computer screen or monitor. In one aspect, the 
invention displays the high frequency data from the QRS 
complex in microvolts adjacent to a display of the conven- 
tional ECG data in millivolts. In another aspect of the inven- 25 
tion, the high frequency data are analyzed for their root mean 
square (RMS) voltage values (as well as for related values 
such as the high frequency energy (HFQE) and high fre- 
quency integral of absolute value (HFAV), both described 
hereinafter), with the discrete, lead-by-lead values being dis- 3Q 
played in real time as useful diagnostic indicators for 
ischemia. Also, the high frequency data from the QRS com- 
plex are analyzed with imbedded algorithms to determine the 
presence or absence of reduced amplitude zones, referred to 
hereinafter as “RAZs”. The given RAZ, of which there are at 35 
least three possible variations (i.e., the “Abboud” RAZ, the 
“NASA” RAZ, and the skewness-kurtosis or “S-K” RAZ, all 
described below), is displayed as a real-time “go, no-go” 
signal on the screen. Finally, in still another aspect of the 
invention, not only the presence or absence any one of the 4Q 
three variations of RAZs, but also the RMS and related values 
(HFQE and HFAV) are displayed against time as time varying 
signals. 

In still another aspect of the invention, the electrocardio- 
graph of this invention detects and aligns R- waves and QRS 45 
complexes and analyzes these ECG signals after digitization 
at high sampling rates of at least 500 samples per second, but 
preferably at sampling rates of 1000 samples per second or 
greater. The system also signal averages consecutive QRS 
electrocardiographic complexes in a user-adjustable fashion 50 
to increase the signal-to-noise ratio. The resulting averaged 
signal is filtered, using non-recursive digital bandpass filters 
with varying high and low frequency cutoffs, and displayed in 
real time along with several other derived numerical measures 
described herein including the power spectrum of the filtered 55 
data. The resulting displays provide a clinician with real-time 
information with respect to changes in the high frequency 
ECG complex that may be indicative of myocardial ischemia, 
myocardial infarction, or of changes in myocardial conduc- 
tion that are unrelated to myocardial ischemia or infarction. 60 

The invention includes a number of features that are neither 
shown nor suggested in the art, including a real-time display 
of cardiac electrical data that have been manipulated in such 
a way as to provide a clear indication of ischemia and infarc- 
tion. The invention further includes a number of user select- 65 
able parameters to enhance the information provided to the 
clinician. Finally, the invention provides a number of dis- 
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plays, juxtaposed in real time, to provide a side-by-side com- 
parison of various aspects of the QRS complex in real time. 

These and other features of the invention will be apparent 
to those of skill in the art from a review of the following 
detailed description along with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic diagram of the overall system of this 
invention. 

FIG. 2 is a schematic diagram of a detail of FIG. 1. 

FIG. 3 is a schematic diagram of the logic carried out by the 
invention. 

FIG. 4 is a real-time screen display, showing characteristic 
data obtained from a healthy subject, including a side-by-side 
display of a standard ECG and a filtered (high frequency) 
ECG. 

FIG. 5 is a real-time screen display, showing characteristic 
data obtained from a patient having known cardiac disease, 
including a side-by-side display of a standard ECG and a 
filtered ECG. 

FIGS. 6 and 7 are real-time screen displays, showing the 
configuration of a QRS detector for a normal, healthy subject 
and for a patient having known cardiac disease, respectively. 

FIGS. 8 and 9 are real-time screen displays of cross corre- 
lation between a running, continuously updated waveform, or 
template, and a sensed waveform to determine departure from 
the template from one heartbeat to the next. 

FIGS. 10 and 11 are real-time screen displays of short term 
trends of various data for a healthy subject and for a patient 
having known cardiac disease, respectively. 

FIGS. 12 and 13 are real-time screen displays of long term 
trends of various data for a healthy subject and for a patient 
having known cardiac disease, respectively 

FIGS. 14 and 15 are real-time screen displays of high 
frequency QRS skewness versus kurtosis for a healthy subject 
and for a patient having known cardiac disease, respectively. 

FIGS. 16 and 17 are real-time screen displays of the spec- 
trum of power for normal and diseased patients, respectively. 

DETAILED DESCRIPTION OF A PREFERRED 
EMBODIMENT 

FIG. 1 shows a simplified, functional, block diagram of a 
real-time high frequency QRS electrocardiograph 10 con- 
structed in accordance with the present invention. The inven- 
tion monitors the cardiac function of a patient with a plurality 
of patient electrodes 12. The electrodes provide measure- 
ments of cardiac electrical function at various contact points 
on the skin of a patient in the conventional manner. For 
example, in the conventional 12 -lead configuration, 10 elec- 
trodes placed upon the skin of the patient in the conventional 
configuration result in eight channels of incoming data. These 
eight channels are in turn translated into 12 leads of data on 
the patient monitor inasmuch as data for one of the bipolar 
limb leads and for all of the augmented unipolar limb leads 
can be derived if data for any two of the bipolar limb leads are 
already known. The analog measurements are coupled to a 
console 14 by way of a communications channel such as for 
example a cable 13. The console components are shown in 
greater detail in FIG. 2. 

The console 14 conditions and digitizes the analog signal 
and provides the digitized signal to a computer 16 by way of 
a communications channel 15, which may preferably be a 
conventional cable or a wireless communication channel by 
radio frequency wave. The structure and function of the com- 
puter is shown and described below in respect of FIG. 3. The 
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computer 16 is programmed to display the ECG signal in real 
time, although the ECG signal may also be stored on a digital 
recording medium 18 over a communications channel 17 for 
later display over a communications channel 17*. 

The computer 16 is coupled to a user interface 20 which 
preferably includes communications devices 22 such as a 
mouse, keyboard, and/or touch screen. The user interface 
further includes a monitor 24 for user controllable graphic 
display of the ECG and various aspects of the signal, a feature 
of the present invention. The computer 16 is coupled to the 
interface 20 by way of bidirectional communications chan- 
nels 19 and 19', for example. The aspects of the graphical 
display are shown in greater detail and described below in 
respect of FIGS. 4 through 17. 

FIG. 2 depicts the structure of the console 14 in greater 
detail. As previously described, the patient is wired with a set 
of electrodes 12, such as for example a conventional set often 
electrodes for a 1 2-lead electrocardiograph to monitor car- 
diac function from different aspects of the patient’ s body. The 
electrodes 12 provide a set of analog electrical signals to the 
console 14, where they are received by a pre-amplifier 30 to 
boost the amplitude of the weak ECG signals. The amplified 
signals are then fed to an anti-aliasing filter 32. The filtered 
signals are then fed to an analog to digital converter 34, where 
the signals are digitized at least the Nyquist rate, preferably 
1 ,000 Hz or greater, to retain all of the information contained 
in the analog signals. The sampled/digitized signals are then 
sent to the computer 16 by an appropriate medium 15. 

The operation of the computer 1 6 is depicted in FIG. 3 . The 
console 14 feeds the digitized ECG signals to the computer 
via a communications channel 15, as previously described. 
The computer 16 also interfaces with the user interface 20. 
The computer 1 6 receives the ECG signals into a device driver 
40, which is simply the interface device and program for the 
console and computer. The device driver 40 provides the ECG 
signals in parallel to an R-wave detection block 42 to syn- 
chronize the system for the start of each heartbeat by locating 
the fiducial point. The following requirements must be met 
for temporal averaging to work effectively. First, the signal of 
interest must be repetitive and relatively invariable. Time 
varying signals, such as ectopic or premature complexes, are 
eliminated before averaging by comparing incoming signals 
against a previously established template through the use of a 
real-time cross -correlating technique. Second, the signal of 
interest must be timelocked to a fiducial point, such as near 
the peak of the QRS complex, that is easily detectable and 
serves as a timing reference for the averaging algorithm. If the 
signal of interest does not have a fixed, temporal relationship 
with the timing reference point, the resultant averaged signal 
will be filtered and distorted due to reference jitter, with 
subsequent loss of the high frequency components. Third, the 
signal of interest and the noise must be independent and 
remain independent during averaging. 

Once the fiducial point has been located for each incoming 
beat, the digitized signals are fed to block 44 where initial 
templates of the QRS complexes are generated. The present 
invention includes running signal averages of the QRS com- 
plexes that constitute templates, with the number of indi- 
vidual beats in the running templates being selectable by the 
user on the user interface 20. The user can also determine the 
percent of cross -correlation between each new incoming beat 
and the templates which will be detected as a departure from 
the norm, and which channel(s) that will be used for the 
cross-correlation functions. 

With the user-selected inputs as just described, the system 
in block 48 runs a real-time cross -correlation function on 
each subsequent incoming beat. In block 50, those beats 
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which are below the threshold set by the user (or set by the 
system default) are rejected, with block 52 accepting only 
well-correlated beats to create the running templates. This 
feature helps to eliminate noisy, unreliable waveforms when 
5 creating the running templates. From this block, the user 
interface displays a continuously updated running total of the 
number of beats accepted and rejected in block 51, a feature 
of the invention. 

Block 54 then aligns the beats by detecting and adjusting 
to for any signal jitter, which may be created by any number of 
well-known factors, including minor inconsistencies in the 
detection of the fiducial point, movement by the patient or 
even respiration by the patient. The aligned, jitter-corrected 
waveform is then fed to block 56, a band-pass filter, prefer- 
15 ably 1 50-250 Hz, to select only the frequencies of interest in 
the waveform. Finally, the band-pass signals are fed to block 
58, where, in the present embodiment, the eight channel to 
1 2-lead conversion is performed. 

From this point in the diagram of FIG. 3, blocks 62,64,68, 
20 70, and 72 describe data which are displayed on the user 
interface 20, as shown in FIGS. 4 through 17. Block 62 shows 
the instantaneous, real-time high frequency filtered QRS sig- 
nal for each lead, updating with each new beat as that beat is 
incorporated into the average beat. Beats that are poorly 
25 cross-correlated are rejected and thus the template averages 
and the display will not be altered by such beats. The display 
is shown by element number 114 in FIG. 4. 

Block 64 determines the unfiltered QRS interval onset and 
offset automatically and in real time. This block receives an 
30 input from block 54 which detected and adjusted for any j itter 
in the well cross-correlated beats. The cross -correlation and 
jitter-correction from block 54 are therefore also displayed on 
the user interface, as shown in FIGS. 8 and 9. Block 68 
describes the display of the instantaneous high frequency 
35 QRS RMS, HFQE and HFAV voltages and power spectrum 
generated for each lead using the unfiltered QRS interval 
onset and offset, and shown in FIG. 4 as element numbers 116 
and 128. 

Block 70 describes the running plots of the RAZs (“go/no- 
40 go”) as well as the voltages for the RMS, HFQE and HFAV of 
the high frequency QRS ECG signal versus time for each lead 
as depicted in FIGS. 10 through 13. FIGS. 10 (healthy sub- 
ject) and 11 (subject with known coronary artery disease) 
depict short-term data trends whereas FIGS. 12 (healthy sub- 
45 ject) and 13 (subject with known coronary artery disease) 
depict longer-term data trends. The horizontal (time) axis 
scale is user-adjustable on both the long and short-term trend 
plots, and can be shown either in beats (as depicted here) or in 
seconds. Clinicians can utilize these trends to assess how a 
50 monitored patient’s cardiac function has changed over time, 
up to and including the present time. Specifically, clinicians 
can identify whether and when RAZs have developed or 
disappeared during the period of monitoring, as well the 
degree to which the RMS and related voltages of the high 
55 frequency QRS complex have changed over the same period 
of monitoring. Realization of such changes is particularly 
valuable to clinicians during situations when the presence or 
absence of incipient myocardial ischemia or infarction needs 
to be immediately identified, when the success or failure of 
60 invasive or noninvasive treatments admini stered for i schemia 
and infarction needs to be immediately recognized, and/or 
when cardiovascular responses during pharmacological or 
exercise stress tests or during a patient’s ambulatory activities 
needs to be assessed. 

65 F inally, block 72 describes another feature of the invention, 

online instantaneous RAZ detection for each lead. The pres- 
ence of a “RAZ”, or “reduced amplitude zone” within the 
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envelope of the averaged high-frequency QRS signal, may be 
an indication of abnormal cardiac function. A RAZ, as origi- 
nally defined by Abboud (but only in the context of off-line 
analyses), classically occurs when at least two local maxima 
of the upper envelope or two local minima of the lower 5 
envelope are present within the high frequency QRS signal. A 
local maximum or minimum is in turn defined as an envelope 
sample point (peak or trough) within the QRS interval 
wherein the absolute value of its voltage exceeds that of the 
three envelope sample peaks immediately preceding and fol- 10 
lowing it. The RAZ is thus the region lying between the two 
neighboring maxima or minima. The present invention per- 
forms a real-time calculation, looking for local maxima and 
minima of the QRS envelope not only according to previously 
published off-line criteria of Abboud (i.e., “RAZ^”, or the 15 
Abboud RAZ) but also separately and especially according to 
new criteria that improve the specificity and accuracy of RAZ 
detection, particularly for the online setting. 

One especially important modification contributing to such 
improvement is the requirement that the second, smaller local 20 
maximum within the high frequency QRS signal that defines 
the presence of the RAZ be at least two times, and optionally 
up to three to four times, larger in amplitude than the RMS 
“noise” of the high frequency signal that is present outside of 
the QRS interval (i . e. , in any given chosen segment within the 25 
PR, ST or TP intervals). Another modification criterion 
implemented optionally by the user into the online RAZ 
detector is that the absolute value of the amplitude of the 
smaller of the two local maxima (or minima) constituting the 
RAZ must be at least a certain percentage of that of the larger 30 
of the two local maxima (or minima ). These modifications as 
optionally implemented by the user within the present inven- 
tion, together with the requirement that at least two local 
maxima of the upper envelope “and” (rather than “or”) two 
local minima of the lower envelope be present to constitute a 35 
RAZ, are collectively referred to as the criteria that define the 
“NASA RAZ”, or “RAZ^”. 

Yet another novel set of user-selectable criteria used in the 
present invention to assess the presence or absence of a “sta- 
tistical” type of RAZ concerns the use of real-time calculation 40 
of the both the skewness and the kurtosis of the incoming high 
frequency QRS signal. When present, the skewness-kurtosis 
type of RAZ is referred to as the RAZ^_^. The presence of a 
RAZ 4 , RAZ^y, and/or RAZ^ by “go/no-go” indicators on 
the display, as shown by element 117 in FIG. 4, may also be 45 
displayed as a running parameter with time, in a manner 
similar to FIGS. 10-13. 

Referring now more particularly to FIGS. 4 and 5, a feature 
of the present invention is the simultaneous display in real- 
time of various aspects of cardiac function with respect to 50 
both the conventional and high frequency QRS ECG. Further, 
the present invention provides simultaneous, side-by-side or 
vertical displays of that data to provide the clinician with a 
tool to compare these aspects of the ECG with one another for 
a more complete picture of cardiac function than has been 55 
previously available. The present invention also displays the 
conventional and high frequency 12-lead configurations on 
one user interface display, while displaying the presence or 
absence of RAZs on the same display to alert the clinician to 
potentially abnormal cardiac function. 60 

The display of FIGS. 4 and 5 includes a conventional, low 
frequency ECG signal, designated by the element number 
112. The display signal 112 includes the signals from leads L, 

II, III, aVR, aVL, aVF, and V1-V6 in the conventional man- 
ner. Further, the display signal 112 slides from left to right, 65 
updated with each beat, also in the conventional manner, as a 
reference for the clinician. Positioned immediately adjacent 
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to or below the conventional ECG display 112 is a display 114 
of a running, instantaneous filtered (i.e. high frequency) QRS 
ECG signal, one for each of the twelve leads, corresponding 
to the individual leads of the display 112. The display signal 
114 includes the signals from leads I, II, III, aVR, aVL, aVF, 
and V1-V6 to correspond to like signals in the display 112. 
Above each high frequency QRS signal is the instantaneous 
RMS value 116 in that particular lead, and an instantaneous 
RAZ indicator 117, defaulted to the RAZ type of the user’s 
choice (in this case the RAZ^). 

At the bottom of the screen display of FIGS. 4 and 5 is a 
tool bar 120. The tool bar 120 provides user control and 
display of various data useful to the clinician. The tool bar 
includes a user selectable indicator 122 in which the user can 
select which of a plurality of QRS complexes to capture. 
Displays 124, 126, and 128 show the timing of the QRS onset, 
QRS offset, and the duration of the QRS complex, respec- 
tively. A display 130 shows the total number of heartbeats 
which have been detected during any particular run, while a 
display 132 shows the number of beats rejected and a display 
134 shows the number of beats whose processing is pending, 
if any. Toggle buttons 136 and 138 permit the user to turn on 
and turn off the cross correlation function and the beat rejec- 
tion function, respectively. Display buttons 140, 142, 144, 
and 146 permit the user to select other displays for the screen, 
as shown. A toggle button 148 permits the user to turn on and 
turn off an autoscale function, and a print button 148 permits 
the user to print a particular screen capture at his discretion. 

Of particular note, in comparing FIGS. 4 and 5, is the small 
number of reduced amplitude zones and lit RAZ indicators 
117 for the healthy patient shown in FIG. 4, versus the much 
larger number of reduced amplitude zones and lit RAZ indi- 
cators for the subject with known coronary artery disease, 
shown in FIG. 5. 

FIGS. 6 and 7 show configuration screen 160 for the online 
QRS interval detector of the invention, selectable by button 
146 of FIGS. 4 and 5. The configuration screen of FIGS. 6 and 
7 permits the user to select a period within the P-R interval of 
the ECG, with specific pre-Q offset and interval width selec- 
tors 162 and 164, respectively, that will aid in the determina- 
tion of the onset of the QRS complex. The values for these 
parameters shown in FIGS. 6 and 7 have been selected based 
on initial experience to provide satisfactory performance of 
the invention for the broadest array of subjects, since cardiac 
function typically varies from patient to patient. Similarly, 
selectors 166 and 168 permit the user to select the offset and 
width, respectively, of a portion of the ST segment that will 
aid in the determination of the offset of the QRS complex. The 
configuration screen 160 also displays various parameters 
measured from the current heartbeat, as shown. An indicator 
170 shows the user which of the leads is selected, in the case 
of FIGS. 6 and 7 that lead is lead II. A RAZ indicator 117 is 
also provided for the convenience of the user. 

FIGS. 6 and 7 also include superimposed displays of a low 
frequency signal averaged ECG 172 and an averaged high 
frequency filtered signal 174. The ordinate on the left of 
FIGS. 6 and 7 shows the averaged low frequency signal 172 in 
millivolts, and the ordinate on the right shows the averaged 
high frequency signal 174 in microvolts. 

FIG. 6 shows a single local maximum 176 and a single 
local minimum 178. However, FIG. 7, depicting the wave- 
forms for a patient with known coronary artery disease, shows 
two local maxima 180 and 182, and two local minima 184 and 
186. Since the maxima and/or minima occur within the QRS 
complex, and are separated by at least three envelope sample 



US 7,539,535 B1 


9 


10 


points with lesser absolute amplitudes between them, they 
define a reduced amplitude zone as also shown on the lit RAZ 
indicator 117. 

FIGS. 8 and 9 depict cross correlation screens 190 of the 
invention. The screen 190 is selected by a user with button 5 
144 (FIG. 4). The screen includes a display 192 of a template, 
averaged over a user selectable number of beats. A display 
194 shows the low frequency ECG of the current beat, and a 
display 196 shows the cross correlation between the wave- 
forms of display 192 and display 194. A user may select the to 
threshold for the cross correlation, below which the beat is 
rejected, by a selector 197 and, if a beat is rejected, an indi- 
cator light 198 illuminates. 

Note that in FIG. 9 the current heartbeat shown on the 
display 194 is noisy and has a different shape than the tern- 15 
plate. The current beat is also shifted in time from the tem- 
plate of display 192. Consequently, the cross correlation 
between them is poor, as reflected in the display 196, which 
only has a peak of 0.95 1 (i.e. below the user selectable thresh- 
old set at 0.970 in this case) and the beat is rejected, as shown 20 
by the indicator light 198. Thus, the current beat is not incor- 
porated into the running template of the averaged ECG. 

FIGS. 10 through 13 depict trend lines for the RMS, 
HFQE, and HFAV, as well as show the current value for these 
parameters. These displays further include lines depicting the 25 
presence versus the absence of three types of RAZ as 
described herein. These displays illustrate another feature of 
the present invention, that of showing a trend of parameters 
overtime. Specifically, the trend lines shown on a screen 200, 
which are continuously updated in real time, illustrate the 30 
presence (versus the absence) of all three types of RAZs over 
a user-selectable period of time (always up to and including 
the present time) as well as the voltage trends for RMS, HFQE 
and HFAV over this same user-selectable period of time for all 
the leads attached to the patient. FIGS. 10 and 11 provide 35 
examples of trend lines on a short-term time scale for a 
healthy subject and for a patient with known coronary artery 
disease, respectively, whereas FIGS. 12 and 13 provide cor- 
responding examples of trend lines on a long-term time scale 
for these same individuals . The short and long-term trend line 40 
plots can be accessed simultaneously by the user, and, to 
reiterate, the time intervals (horizontal axes depicted here in 
units of beats) on both plots are completely selectable by the 
user. The particular parameters depicted on these plots are 
intended to be illustrative only, and other parameters may be 45 
similarly illustrated. The ordinate for each plot is in micro- 
volts with respect to the RMS, BFQE and HFAV (the RAZs 
are unitless, “go versus no-go” entities), whereas the abscissa 
can be in either beats or seconds, for characterization of the 
trend. 50 

Numerical measures of the high frequency QRS ECG may 
be calculated in several ways. These measures are important 
because they often decrease when ischemia is present. Per- 
haps the most popular measure is the root mean square (RMS) 
voltage of the QRS signal, which is equivalent to the “area 55 
under the curve” of the power spectrum, defined as 


RMS = 


Z V 

A i=ufqon 

\ UFQRSD 


60 


where X z is the filtered voltage at a given sampling point, 65 
ufqon and ufqoff are the onset and offset, respectively, of the 
QRS interval, and UFQRSD is the unfiltered QRS interval 


duration as defined by ufqon and ufqoff. In this context, the 
term “onset” means the start of the QRS interval, and the term 
“offset” means the end of the QRS interval. This is the pri- 
mary numerical measure used in preferred embodiment of the 
present invention. 

Other numerical measures of the high frequency QRS sig- 
nal have also been proposed (and used strictly in an off-line 
fashion) by Xue et al. (see Xue, Q., B. R. Reddy, and T. 
Aversano. Analysis of high- frequency signal-averaged ECG 
measurements. J Electrocardiol 28: 23945, 1995). These 
numerical measures include the high frequency integral of 
absolute values (HFAV) and the high frequency QRS energy 
(HFQE). Xue et al. have defined HFAV and HFQE as follows: 


uqoff+iOms 

HFAV = ^ \Xi-AVNL\ 

i=uqon-l0ms 


and, 


uqoff +I0ms 

HFQE = Yj ( Xi-AVNL f 

i=uqon—l0ms 


wherein AVNL equals the average noise level of the filtered 
signal in the ST segment in a 40 ms window located 60 ms 
from the QRS offset. It should be noted that in their own 
definitions for HFAV and HFQE (as shown above), Xue et al. 
“pad” both the QRS onset and offset by extra 1 0 ms each in an 
effort to reduce noise and noise variability, presumably to 
compensate for potential inaccuracies and inconsistencies 
related to the determination of the QRS interval. 

In the presently preferred embodiment of the present 
invention, the definitions for HFAV and HFQE are modified 
in two ways. First, because the present invention provides a 
means of viewing the high and low frequency ECG signal in 
real time, thereby providing reliability for the determination 
of the QRS interval, the necessity (or lack thereof) of using 
the 10 millisecond padding periods is left to the discretion of 
the user. Second, the present invention preferably uses the PR 
interval or TP segment rather than the ST segment to deter- 
mine the AVNL of the baseline, since a segment of the cardiac 
cycle wherein neither depolarization nor repolarization is 
present is preferred. In the device of the present invention 
disclosed herein, AVNL is determined as the RMS noise of 
the filtered signal that is within a 25 ms interval in the PR 
segment. 

In the preferred embodiment, HFAV is measured as: 


uqoff 

HFAV = Yj \Xi-AVNL\ 

\-uqon 


wherein AVNL (average noise level) equals the “noise” of the 
high frequency (i.e., filtered) signal within that portion of the 
PR interval just noted. 

Further, high frequency QRS energy (HFQE) is calculated 
as: 
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uqoff 

HFQE = ( Xi-AVNL ) 2 

i=uqon 


wherein AVNL is determined in the same way as for HFAV. 

A useful characterization of a set of data includes skewness 
and kurtosis. Skewness is a measure of symmetry, or more 10 
accurately, the lack of symmetry. A distribution, or data set, is 
symmetric if it looks the same to the left and right of the center 
point. Kurtosis is a measure of whether the data are peaked or 
flat relative to a normal distribution. That is, data sets with a 
high kurtosis tend to have a distinct peak near the mean, 15 
decline rather rapidly, and have heavy tails. Data sets with low 
kurtosis tend to have a flat top near the mean rather than a 
sharp peak. Auniform distributionwouldbethe extreme case. 

By displaying skewness values alone, kurtosis values alone, 
or skewness versus kurtosis graphically with one value on the 20 
ordinate and the other on the abscissa, the clinician may 
obtain information from the high frequency QRS signals 
potentially indicative of cardiac disease. Referring to FIGS. 

14 and 15, a display 250 is provided to show real-time, con- 
tinuously updated plots of skewness versus kurtosis against a 25 
background line of a normal population distribution function. 

If a patient’s S-K data in any given lead fall on or above that 
lead’ s distribution function, the data for that lead are consid- 
ered “normal”, as shown in most of the leads of FIG. 14, 
whereas data from any given lead falling below that lead’s 30 
distribution function, as shown in most of the leads of FIG. 

15, provides a visual indication of potential cardiac disease, 

i.e., a positive “RAZ^.^”. 

Referring now to FIGS. 16 and 17, a display 270 is pro- 
vided to show the power spectrum of the high frequency QRS 35 
signal for each lead. A change in the shape of the spectral plot 
such that two distinct peaks occur rather than one may be 
further indication of cardiac malfunction. 

The principles, preferred embodiment, and mode of opera- 
tion of the present invention have been described in the fore- 40 
going specification. This invention is not to be construed as 
limited to the particular forms disclosed, since these are 
regarded as illustrative rather than restrictive. Moreover, 
variations and changes may be made by those skilled in the art 
without departing from the spirit of the invention. 45 

We claim: 

1. A process of operating a computer, comprised of a gen- 
eral purpose data processor of known type and a display 
device, to enable the data processor to execute a plurality of 50 
steps in an object program, such that the same results will be 
produced when using the same given data, comprising the 
steps of: 

a. inputting a patient’s electrocardiograph signal to the 
computer via a communications channel, the signal 55 
comprising a succession of waves; 

b. using the data processor to independently isolate the 
QRS interval from each of the waves of the signal; 

c. using the data processor to independently analyze data 

within each of the isolated 60 

QRS intervals relating to cardiac function of a patient in 
real-time, the step of analyzing including band-pass fil- 
tering the data between about 150 Hz and about 250 Hz; 
and 

d. using the display device to display the analyzed data, the 65 
signal, the waves of the signal, the QRS intervals, or any 
combination beat -by-beat in real-time. 
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2. The process of claim 1, wherein the step of analyzing the 
data within the QRS interval further comprises the steps of: 

using the data processor to generate an initial template of 
the signal using the first QRS interval from the first wave 
of the signal; 

using the data processor to cross-correlate each subsequent 
independently isolated QRS interval against predeter- 
mined criteria wherein accepted QRS intervals are 
incorporated into the initial template to create a running 
template in real-time; 

using the data processor to align each independently iso- 
lated QRS interval with the running template in real- 
time; and 

using the data processor to calculate predetermined, select- 
able characteristics of each QRS interval in a beat -by- 
beat, real-time manner. 

3 . The process of claim 2, wherein the step of using the data 
processor to calculate is comprised of calculating a reduced 
amplitude zone, 

wherein the reduced amplitude zone occurs when at least 
two local maxima or at least two local minima are 
present in the filtered data; and 

wherein each local maximum or each local minimum is a 
peak or a trough, respectively and wherein the absolute 
value of each peak or trough’s voltage exceeds that of a 
predetermined number of sample points within the fil- 
tered data immediately preceding and following each 
peak or trough. 

4. The process of claim 3, wherein the absolute value of the 
amplitude of the smallest of the at least two local maxima or 
the smallest of the at least two minima is at least a predeter- 
mined percentage relative to the largest of the at least two 
maxima or the largest of the at least two minima, respectively. 

5. The method of claim 3, wherein the absolute value of the 
amplitude of the smallest of the at least two local maxima or 
the smallest of the at least two local minima is at least a 
predetermined multiple of the high frequency noise level of a 
predetermined iso-electric portion. 

6 . The method of claim 2, wherein the step of using the data 
processor to calculate is comprised of calculating a skewness, 
kurtosis, or both of the filtered data. 

7. The process of claim 1, wherein the step of using the data 
processor to independently analyze is further comprised of 
using the data processor to compute the presence or absence 
of reduced amplitude zones from the high frequency filtered 
data associated with the QRS interval. 

8. A process of operating a computer, comprised of a gen- 
eral purpose data processor of known type, to enable the data 
processor to execute a plurality of steps in an object program, 
such that the same results will be produced when using the 
same given data, comprising the steps of: 

a. inputting a patient’s electrocardiograph signal to the 
computer via a communications channel, the signal 
comprising a succession of waves; 

b. using the data processor to independently isolate the 
QRS interval from each of the waves of the signal; 

c. using the data processor to independently analyze data 
within each of the QRS intervals relating to cardiac 
function of a patient in real-time, the step of analyzing 
including band-pass filtering the data between about 1 50 
Hz and about 250 Hz; 

d. connecting a storage device to the data processor; and 

e. using the storage device to store the analyzed data, the 
signal, the waves of the signal, the QRS intervals, or any 
combination. 

9. An electronic storage database of the analyzed data, the 
signal, the waves of the signal, the QRS intervals, or any 
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combination resulting from the process described in claim 8 * * * * * 

wherein the database is stored in the storage device. 



